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The inclusive decay B —> X s y is well known as one of the most important flavour observables 
within the indirect search for new physics [jl]]. The present experimental accuracy already reached 
the 10% level, as reflected in the world average of the present measurements [2]: 

BR[5 -> X s y] = (3.39") x 1(T 4 . (1) 

In the near future, more precise data on this mode are expected from the /^-factories. Thus, it 
is mandatory to reduce the present theoretical uncertainty accordingly. Non-perturbative effects 
are naturally small within inclusive modes [[!]] ; also additional non-perturbative corrections due to 
necessary cuts in the photon energy spectrum are under control (see [||]). As was first noticed in 
[4j], there exists a large uncertainty in the theoretical NLL prediction related to the renormalization 
scheme of the charm-quark mass on which we focus in this article. The reason is that the matrix 
elements (sy\0\p\b), through which the charm-quark mass dependence dominantly enters, vanish 
at the lowest order (LL) and, as a consequence, the charm-quark mass does not get renormalized in 
a NLL calculation, which means that the symbol m c can be identified with m C)PO i e or with the MS 
mass m c (p c ) at some scale p c or with some other definition of m c . In a recent theoretical update 
of the NLL prediction of this branching fraction the ratio m c /mb was varied in the conservative 
range 0.18 < m c /mi, < 0.31 that covers both the pole mass value (with its numerical error) and the 
running mass m c (p c ) value (with p c £ [m c ,m b ]), leading to [§]: 

BR[B -► X s y\ = (3.70±0.35| mfM ±0.02|cKM±0.25| param . ±0.15| sca i e ) x 1CT 4 . (2) 

The only way to resolve this scheme ambiguity in a satisfactory way is to perform a systematic 
NNLL calculation. Working to next-to-next-to-leading-log (NNLL) precision means that one is 
resumming all the terms of the form 

(a s (m b )) p a n s (m b )\og n (m b /M), (p = 0,1,2). (3) 

where M = m t or M = mw, n = 0,1,2,... . Such a calculation is most suitably done in the frame¬ 
work of an effective low-energy theory. The effective interaction Hamiltonian can be written as 

H efi = -4G F /V2 V tb V t * ^Ci(p.M) Oi{p), (4) 

where Oj(p) are the relevant dimension 6 operators and Cj(p.M) are the Wilson coefficients. 

Parts of the three principal calculational steps leading to the NNLL result within the effective 
field theory approach are already done: (a) The full SM theory has to be matched with the effective 
theory at the scale p = pw, where pw denotes a scale of order mw or m t . The Wilson coefficients 
Ci(pw ) only pick up small QCD corrections, which can be calculated in fixed-order perturbation 
theory. In the NNLL program, the matching has to be worked out at the order oq. The matching 
calculation to this precision is already finished, including the most difficult piece, the three-loop 
matching of the operators Ojfi [6j. (b) The evolution of these Wilson coefficients from p = pw 
down to p = pi, then has to be performed with the help of the renormalization group, where p b 
is of the order of m b . As the matrix elements of the operators evaluated at the low scale p b arc 
free of large logarithms, the latter are contained in resummed form in the Wilson coefficients. For 
the NNLL calculation, this RGE step has to be done using the anomalous-dimension matrix up to 
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order a/. While the three-loop mixing among the four-quark operators O, (/' = 1.... . 6 ) [7j] and 
among the dipole operators Oq% [ Hi] are already available, the four-loop mixing of the four-quark 
into the dipole operators is still an open issue, (c) To achieve NNLL precision, the matrix elements 
(XsY| Oi(pb)\b) have to be calculated to order a 2 precision. This includes also bremsstrahlung 
corrections. In 2003, the (a 2 rif) corrections to the matrix elements of the operators 0\,0i,0-i,0% 
were calculated [|9 ]. Complete order a 2 s results are available to the ( Oq . O-j) contribution to the 
decay width [jH]]. Recently, also order a 2 terms to the photon energy spectrum (away from the 
endpoint £.™ ax ) were worked out for the operator (9 7 [|TTjl. 

In ref. [12] a strong motivation for this complicated NNLL effort was given by calculat¬ 
ing those NNLL terms that are induced by renormalizing the charm-quark mass in the NLL ex¬ 
pressions, i.e. those terms that are sensitive to the definition of the charm-quark mass. These 
terms correspond to 8m c insertions in the diagrams related to the NLL order matrix elements 
M \ l 2 (m c ) = (sy\O ia (Bb)\b) and Mj r | ms (m c ) = (syg\Oi 2 (pb)\b) (for an example, see the left di¬ 
agram in Fig. []]). 




Figure 1: Left: Typical 8 m c insertion diagram. Right: Typical diagram with a self energy insertion. 

The sum 8 M™ 1 ^ (m c ) ■ 8m c of all these insertions can be obtained by replacing m c by m c + 8m c 
in the 0 (a|) results M™^ E \m c ), followed by expanding in 8m c up to linear order: 

M\f z) ( m c + 8 m c ) = < 2 rt(e) (, m c ) + ( m c ) ■ 8m c + 0((8m c ) 2 ). (5) 

As 8m c is ultraviolet-divergent, the matrix elements (m c ) are needed in our application up to 

order s 1 , as indicated by the notation in eq. (!5i). In Jl2| ] the explicit analytical results for these matrix 
elements are given in such a way that they can be used in a future complete NNLL calculation. The 
explicit shift 8m c depends of course on the renormalization scheme. When aiming at expressing 
the results for (, m c ) in terms of m, (ph) or m c>po i e , the shift reads (C 7 = 4/3) 

Om c yfUjy J — ^ ^ ^^c,pole — Cp ( ~ ~\~ 3 111 —2 4 ) Jfl c ,pole • 

The infinities induced by the 1/e terms in 8m, get cancelled in a full NNLL calculation, in particular 
by self-energy diagrams depicted in the right diagram in Fig. [j]. When implementing these self¬ 
energy insertions, we only took into account the Z] (p 2 = m 2 ) piece, i.e. that part of the one- 
loop self-energy which only gets renormalized by the mass parameter. When used at the fixed 
momentum p 2 = m 2 , this piece is gauge-independent. 
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Our final estimates are given in Fig. ||for three different values of p/„ where pi, represents the 
usual renormalization scale of the effective field theory. Within each vertical string, the solid dot 
represents the branching ratio using the pole mass m c>p0 i e , while the open symbols correspond to 
the MS mass m c {p c ) for Be = 1-25 GeV (triangle), p c = 2.5 GeV (quadrangle) and p c = 5.0 GeV 
(pentagon). For each pi, the left string shows the value of the branching ratio at the NLL level, 
while the right string shows the corresponding value where, in addition. 8 m c mass insertions and 
X ] (p 2 = ml) insertions were taken into account. Because the combination of these insertions is 
zero by construction for the pole scheme, the solid dots are at the same place in the left and the 
right string for a given value of p/,. We stress that all the statements made in the following are 
independent of this absolute normalization introduced by the additional £i {p 2 = ml) insertions, 
because we refer to the reduction of the error only. From Fig. [2] we see that the error related to 
the charm-quark mass definition is significantly reduced when the NNLL terms connected with 
mass insertions are taken into account. Taking as an example the results for pb = 5 GeV, we find 
that at the NLL level the branching ratio evaluated for m c ( 2.5 GeV) is 12.6% higher than the one 
based on m CiPO i e . Including the new contributions, these 12.6% get reduced to 5.1%. One also 
can read off an analogous significant reduction within the MS scheme itself. However, to obtain a 
NNLL prediction for the central value of the branching ratio, it is of course necessary to calculate 
all NNLL terms. 



M b [GeV] 


Figure 2: BR(i? —> X s y) for three values of pi, (see text for more details). 
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